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1. INTRODUCTION 

Bose-Einstein condensation in atomic hydrogen was observed for the first time just a 
few weeks before this session of the Enrico Fermi School, and so these lectures constitute 
a somewhat breathless first report. However, the search for Bose-Einstein condensation 
(BEC) in hydrogen began many years ago, and it has a long and colorful history. Out 
of that history emerged an extensive body of knowledge on the behavior of hydrogen at 
low temperatures that provided the foundation for achieving BEC in hydrogen. Much 
of it has been described in reviews [1-4]. Consequently, we will dwell on only those 
features of that research that provide essential background, and concentrate on the most 
recent developments in which BEC in hydrogen advanced from a long period of being 
tantalizingly close, to being real. 

2. ORIGINS OF THE SEARCH FOR BEC IN AN ATOMIC GAS 

In April, 1976, W. C. Stwalley and L. H. Nosanow published a letter summarizing 
studies on the equation of state of spin-polarized hydrogen, H{ [5]. Because there are 
no bound states of molecular hydrogen in the triplet state, H{ behaves like a simple 
monatomic gas, but a gas with a remarkable property. Because of the weak H[ - H{ 
potential and the atom's low mass, H[ remains a gas at temperatures down to T = 0. 
Consequently, it might be possible to cool H{ to the quantum regime and achieve BEC. 
That paper essentially launched the search for BEC in an atomic gas. It triggered a 
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flurry of experiments and enough activity for an entire session of the December, 1978, 
APS meeting to be devoted to the stabilization of hydrogen [6]. 

The critical density for the BEC transition in a non-interacting gas is n c = 2.612A^ 3 , 

where At = \J 2nK 2 /feTm is the thermal de Broglie wavelength. Because n c depends 
on the product of temperature and mass, for a given density hydrogen condenses at a 
higher temperature than any other atom. Hydrogen offered two other attractions as a 
candidate for BEC. The hydrogen atom is generally appealing for basic studies because 
its structure and interactions can be calculated from first principles. Furthermore, H\, 
constitutes a nearly ideal Bose gas: it has an anomalously small s-wave scattering length 
and its interactions are weak. 

In the years since the search for BEC in atoms started, there was a revolution in 
techniques for cooling and manipulating atoms using laser-based methods that culmi- 
nated in the achievement of BEC in alkali metal atoms [7-9] . In light of these advances, 
hydrogen's special attractions must be viewed from a new perspective. Spin-polarized 
hydrogen's unique property of remaining a gas at T — is evidently not essential for 
BEC. Although all species except helium should be solid at sub-kelvin temperatures, 
laser-cooled atoms do not even liquefy because they never hit surfaces. In the absence 
of surface collisions, a gas to liquid transition requires three-body collisions to initiate 
nuclcation. At the densities used to achieve BEC in alkali metal atom systems, however, 
the three-body recombination rate is so low that it can be neglected at all but the high- 
est densities. Hydrogen's relatively high condensation temperature is also not a crucial 
advantage. Laser cooling methods make it is possible to cool alkali metal gases to the 
microkelvin regime, far colder than possible by conventional cryogenic means. Once the 
atoms have achieved the laser-cooling temperature limit, they can be efficiently cooled 
into the nanokelvin regime by evaporative cooling. Finally, hydrogen's close to ideal 
behavior must now be regarded as a serious experimental disadvantage. This is because 
all routes to BEC used so far involve evaporative cooling. Evaporation requires collisions 
for maintaining thermal equilibrium as the system cools. Because of hydrogen's small 
scattering length, its collision cross section is tiny and evaporation is much slower than 
in other systems. 

Nevertheless, now that hydrogen can be Bose-Einstein condensed, its simplicity con- 
tinues to give the atom unique interest. The techniques, condensate size, and the general 
conditions for condensation are different from those of other realizations and one can 
expect that this new condensate will open complementary lines of research. 

Spin polarized hydrogen is created by the magnetic state selection of hydrogen at 
cryogenic temperatures. In high magnetic fields the electron and proton spin quantum 
numbers are m e — — |, m p = ±| for HI and m e = +|, m p = ±5 for H]. The 
governing parameter in magnetic state selection is the interaction energy in a strong 
magnetic field, B. In temperature units, this is T a = fi B B/k B = 0.67(£?/tesla) K. Here 

is the Bohr magneton, and fee is Boltzmann's constant. For a field of 10 T, which is 
readily achieved in a superconducting solenoid, T = 6.7 K. At a temperature of 0.3 K 
the ratio of densities n{H]) / n(i?|) ~ exp(-2T /T), which is about 10~ 20 . The spin 
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polarization is essentially 100%. 

Another crucial experimental parameter in creating H[ is the binding energy for ad- 
sorption on a liquid helium surface, E\>. Hydrogen must make collisions with a cold 
surface in order to thermalize at cryogenic temperatures, but if the atoms become ad- 
sorbed the gas will rapidly recombine. In thermal equilibrium, the surface density a 
and volume density n of H[ are related by a — At« exp (— E\> / k^T) . For hydrogen- 
helium, Eb = 1 K. Consequently, for temperatures below about 0.1 K the atoms move 
to the walls where they can recombine by two- and three-body processes. The initial 
searches for BEC were in the temperature regime 0.2 <~ 0.7 K. The critical density at a 
temperature of 0.5 K is n c = 6.8 x 10 19 cm~ 3 . 

Spin-polarized hydrogen was first stabilized by I. F. Silvera and J. T. M. Walraven in 
1980 [10], and magnetically confined by the M.I.T. group [11]. In these experiments the 
field of a superconducting solenoid provided both state-selection and axial confinement. 
Radial confinement was provided by a supcrfluid helium-coated surface. 

The highest density achieved with HI under controlled conditions was 4.5 x 10 18 cm~ 3 , 
at a temperature of 0.55 K [3]. The density was limited by three-body recombination, 
a process that had been predicted by Kagan et al. [12]. Because the heat generation in 
three-body recombination increases as the cube of the density, at higher density the gas 
would essentially self-destruct. An alternative route to BEC was required. 

3. HYDROGEN TRAPPING AND COOLING 

The new route toward BEC led to temperatures in the microkelvin regime where n c 
is so low that three-body recombination is unimportant. At temperatures much below 
0.1 K, however, surface adsorption and recombination become prohibitive. To avoid 
surfaces, Hess suggested strategies for confining HI atoms (in the "low field seeking" 
states) in a magnetic trap, and cooling the gas by evaporation [13]. 

Trapping and cooling requires an irreversible process for losing energy. In laser- 
cooling and trapping experiments, this process is spontaneous emission. Unfortunately, 
laser methods are not well suited to hydrogen. Aside from the lack of convenient light 
sources, the laser cooling temperature limit for hydrogen is relatively high. The limit is 
determined by the recoil energy for single photon emission, and in hydrogen it is more 
than a millikelvin. The method proposed by Hess employed elastic scattering for both 
trapping and cooling. 

The magnetic trap proposed by Hess consists of a long quadrupolc field to confine the 
atoms radially with axial solenoids at each end to provide axial confinement — an elon- 
gated variant of the "Ioffc-Pritchard" configuration [14]. The Ioffe-Pritchard potential 
is 

(i) v(r) = V(M 2 + (/^ 2 + e? - e, 

with radial potential gradient a, axial curvature 2/3, and bias energy 9. At low energies 
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the potential is quasi-harmonic with radial and axial oscillation frequencies 



(2a) 



a 



Vm^z 2 + 9) 



(26) 




The trapping fields are produced inside a cell that confines the gas cloud while atoms 
are loaded into the trap. The walls of the helium-coated cell are held at a temperature 
of about 275 mK. The cell is filled with a puff of hydrogen and helium from a low 
temperature discharge and the walls of the cell are then cooled, making them "sticky." 
Atoms with high energy leave the trap and stick to the walls. Thermal contact between 
the walls and the trapped gas is broken because these atoms are unable to return. The 
trapped gas cools by evaporation and reaches an equilibrium temperature of about l/12th 
the trap depth. 

To further cool the gas, Hess proposed a process of forced evaporation. The height of 
the potential barrier is slowly decreased, allowing energetic atoms to escape and thus re- 
ducing the average energy of the system. As the system's reduced energy is redistributed 
by collisions, the temperature falls. Unlike ordinary evaporation which has a tempera- 
ture limit determined by the vapor pressure of the material, forced evaporation can be 
continued to almost arbitrarily low temperatures. The process is surprisingly efficient 
because the escaping atoms carry away a great deal of energy, (5 <~ IO^bT per atom. 

Spin-polarized hydrogen was confined by a pure magnetic trap in f 987 by the M.l.T. 
group [15] and also in Amsterdam [16]. Shortly thereafter, in the first demonstration of 
forced evaporative cooling the M.l.T. group achieved a temperature of 3 mK [17]. 

In all of these experiments, the atoms were studied by monitoring the hydrogen flux 
as the atoms were dumped from the trap. With this technique [17,18], the field of one of 
the axial confining coils is reduced, allowing the atoms to escape from the trap. Once out 
of the trap, the hydrogen rapidly adsorbs on the walls of the cell and recombines. This 
heat of molecular recombination — 4.6 eV per event — is measured by a small bolometer 
within the cell but outside the trap. If the confinement field is reduced rapidly compared 
to the thermalization time, then the atom flux as a function of barrier height reveals the 
energy distribution of the gas. From this the temperature can be inferred. A typical 
energy distribution is shown in fig. 1. Integrating the signal gives a measure of the total 
number N of atoms trapped. 

The density can also be determined by observing the decay of the trapped gas. The 
primary decay mechanism in H\ is dipolar relaxation, a process in which the spin angular 
momentum of a pair of atoms is transferred to their orbital angular momentum, while 
one or both of the atoms makes a transition to an untrapped state and escapes. Because 
dipolar relaxation is a two-body process, the density decays according to h = —gn 2 , where 
the dipolar decay constant g has been calculated [19,20] and also measured [15,16]. The 
calculated value is g = 1.2 x 10~ 15 cm 3 s _1 [21]. The total sample decays according to 
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barrier height [mK] 

Fig. 1. - Energy distribution of trapped H] after evaporative cooling: plot of the atom flux 
as the trap barrier is rapidly lowered. The atom flux is measured by detecting the molecular 
recombination heat on a sensitive bolometer. Calculated distributions for three temperatures 
are shown. 



N = —KgN where 

_ fexp{-2V(r)/k B T}dV 
{> Jcxp[-V(r)/k B T]dV 

For dipolar decay the number of trapped atoms decays according to 
(4) N(0)/N(t) = l + ngn t. 

Hence, the initial peak density n D can be extracted from a plot of N(0)/N(t). An example 
of such a plot is shown in fig. 2. Typically k <~ 1/5 and the characteristic decay time at 
a density of 10 14 cm~ 3 is 40 s. 

Dipolar relaxation takes place predominantly where the density is high. This is at 
the minimum of the trap where the mean energy is low. Just as evaporation cools 
by removing the most energetic atoms, dipolar relaxation heats by removing the least 
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Fig. 2. - Determination of the sample density by observing decay due to dipolar relaxation. 
Five samples were indentically prepared and then held for varying times (horizontal axis) before 
being dumped from the trap. The integrated recombination signal from the dump gives the 
number of atoms in the trap, and is plotted on the vertical axis. The straight line fit indicates 
a density of 6.0 x 10 13 cm -3 . From [50] 



energetic atoms. Consequently, the trapped gas comes to a thermal equilibrium in which 
cooling due to evaporation is balanced by heating due to dipolar relaxation. 

These attempts to reach BEC in H] stopped short of the quantum degenerate regime 
by about a factor of six in phase space density [22]. Several complications arose. The 
diagnostic technique of dumping the atoms out of the trap began to give ambiguous 
results because the energy distributions of the samples would change significantly during 
the dumping process. The sensitivity of the detection process was not sufficient. The 
efficiency of the evaporative cooling process degraded significantly, precluding further 
cooling. Finally, even if the quantum degnerate regime could be reached, a condensate 
would degenerate during the sample dumping process, and would have been unobservable. 
These complications required a new diagnostic technique and a more efficient evaporation 
process. 
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4. OPTICAL DETECTION OF TRAPPED HYDROGEN 

Trapped hydrogen can be detected in situ by photoabsorption using one- and two- 
photon transitions. The Amsterdam group has used both methods. Using a Lyman-a 
light source they observed absorption of the principal transition [23] . However, the ab- 
sorption spectrum has a large natural linewidth and displays structure due to the inho- 
mogeneous magnetic field, limiting its use for analyzing momentum. They also observed 
the 15-35 and 15-31? transitions using resonantly enhanced two-photon excitation ex- 
ploiting a virtual state near 2P [24]. The 15-35 transition is insensitive to magnetic 
fields and is potentially well suited for analyzing momcmtum at temperatures down to 
1 /xK. The 15-3D transition has a matrix element that is ten times larger, but a natural 
linewidth that is ten times broader. 

We have employed two-photon Dopplcr-frcc spectroscopy of the 15-25 transition, us- 
ing two-photons from a single laser tuned to 243 nm, twice the Lyman-a wavelength. 
The spectrum is essentially unperturbed by the magnetic field. This method provides 
excellent momentum resolution but suffers from the low excitation rate of a forbidden 
transition. In two-photon spectroscopy using a single light source normally only the 
narrow Doppler-free spectrum is observed, excited by counter-propagating laser beams 
that eliminate the first order Dopplcr effect. However, at very low temperatures in our 
experiment the broad Doppler-sensitive spectrum is also visible, excited by absorbing two 
photons from the same laser beam. Because the momentum distribution in a condensate 
is much narrower than in a normal gas at the same temperature, the Doppler-sensitive 
spectrum is well suited to probing the condensate. In addition to its application to 
observing BEC, Doppler-free spectroscopy of cold trapped hydrogen has a potential ap- 
plication to the precision spectroscopy of hydrogen, in which the 15-25 transition plays 
a central role [25]. Under suitable conditions, time-of-flight broadening disappears and 
a spectral resolution close to the natural linewidth of 1.3 Hz should be possible [26]. 

In our method [27], the atoms are excited by a pulse of 243 nm radiation. For 
observing BEC the pulse length is typically 0.4 ms, but for spectroscopic studies pulses 
as long as 5 ms have been used. (In the absence of an electric field we have observed 
that the 25 atoms live for close to the natural lifetime, 122 ms [27].) Because of the low 
excitation rate, it is not feasible to observe photoabsorption. Instead, the excited atoms 
are detected by switching on an electric field which Stark-quenches the 25 state, mixing 
it with the 2P state which promptly decays. The emitted Lyman-a photon is detected. 
A diagram of the apparatus is shown in fig. 3. The energy equation for two-photon 
excitation of an atom from state i to /, with initial momentum and final momentum 
Pf — Pi + ^(ki + k 2 ), where ki and k 2 are the wave vectors of the laser beams, is 



where the rest mass of the atom in the initial state is m and v a is the unperturbed 
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Fig. 3. - Schematic diagram of the apparatus. The superconducting magnetic coils create 
trapping potential that confines atoms near the focus of the 243 nm laser beam. The beam 
is focused to a 50 /im waist radius and retroreflected to produce the standing wave required 
for Doppler-free two-photon absorption. After excitation, fluorescence is induced by an applied 
electric field. A small fraction of the 122 nm fluorescence photons are counted on a microchannel 
plate detector. The bolometer is used for diagnostic purposes. Not shown is the trapping cell 
which surrounds the sample and is thermally anchored to a dilution refrigerator. The actual 
trap is longer and narrower than indicated in the diagram. 



transition frequency. Expanding, we obtain 

(ki+k 2 )- Pi fi(ki + k 2 ) 2 

(6) V = V ° + . 4nm -t 1 - 6 ^— 8^~ (1 - £) , 

S v ' S v ' 

Av m Au R 

Here Avdi and Avui are the first and second order Doppler shifts, respectively, Av>r 
is the recoil shift, and e = 2hv /mc 2 = 1.1 x 10~ 8 is a relativistic correction which 
accounts for the mass change of the atom upon absorbing energy 2hv . For hydrogen in 
the submillikelvin regime, Ai> D2 <C 1 Hz, and can be neglected. In the Doppler-sensitive 
configuration, ki = k 2 and Az/r = 6.7 MHz. (All frequencies are referenced to the 
243 nm laser source). At a temperature of 50 /uK, Avdi ~ 2.6 MHz. In the Doppler- 
free configuration, ki = — k 2 , and there is no recoil or first order Doppler broadening. 
Doppler-free excitation is achieved by retro-reflecting the laser beam. Nevertheless, in 
this configuration the Doppler-sensitive line is also excited, with the atom absorbing two 
photons from a single laser beam. 



VoPt 




+o 



hv 



BOSE-EINSTEIN CONDENSATION OF ATOMIC HYDROGEN 9 




-60 -40 -20 20 40 60 
laser detuning [kHz at 243 nm] 

Fig. 4. - Spectra of low density samples at a series of temperatures. The linewidth is given 
by the finite interaction time as an atom traverses the Gaussian laser beam. The exponential 
lineshape arises from an average over a Maxwell-Boltzmann distribution of velocities. As the 
sample cools, the characteristic velocity decreases and the line narrows ~ VT. The laser power 
was about 7 mW, and densities were in the range 10 12 ~ 10 13 cm . 



Far from quantum degeneracy, the lineshape for Doppler-scnsitive excitation is the 
familiar Gaussian curve characteristic of a Maxwell-Boltzmann distribution. The shape 
for Doppler-free excitation, however, is quite different — a cusp-shaped double exponential 
[28]: I(y — v ) <~ exp(— \ v — v Q \ /5 ). The linewidth parameter S Q is determined by the 
time of flight of an atom across the laser beam: 6 = u/2nd where u = \j2k^TJrn is the 
most probable velocity and d a is the waist diameter of the laser beam. This expression 
is valid for an untrapped gas far from quantum degeneracy. It neglects collisions and 
other broadening mechanisms, and the natural linewidth of 1.3 Hz. An example of this 
Dopplcr-frcc lineshape is shown in fig. 4. A panoramic spectrum showing the Doppler-free 
and the recoil-shifted Doppler-sensitive lines is shown in fig. 5. 

If the atoms are confined in a radially harmonic trap, the Doppler-free spectrum 
consists of a central line at frequency v a plus a series of sidebands spaced by twice the 
trap frequency, lying under the exponential curve [29] . The intensity of the sidebands is 
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Fig. 5. - Composite 1S-2S two-photon spectrum of trapped hydrogen before condensation. The 
intense, narrow peak arises from absorption of counter-propagating photons by the normal gas, 
and exhibits no first-order Doppler broadening. The wide, low feature on the right is from 
absorption of co-propagating photons. The solid line is the recoil-shifted, Doppler-broadened, 
Gaussian lineshape of the normal gas corresponding to T — 40 fiK. Zero detuning is taken for 
unperturbed atoms excited Doppler-free. All frequencies are referenced to the 243 nm excitation 
radiation. 



governed by the ratio of the atom cloud diameter to laser beam diameter. If this ratio is 
less than 1, only the cental line is excited, an example of Dicke narrowing. 

5. EVAPORATIVE COOLING 

The physical principles and experimental considerations of evaporative cooling have 
been described in detail by Ketterle and van Druten [30]. We summarize here some of 
the principal points. 

Evaporative cooling occurs when highly energetic atoms are permitted to escape from 
a trap at a rate that is kept sufficiently low for the remaining gas to maintain a quasi- 
thermal equilibrium. In this situation the energy distribution is thermal except that it is 
truncated at the depth of the trap, Vt rap . The crucial parameters governing evaporative 
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species 


m (amu) 


a (nm) 


CT e i (cm 2 ) 


T c (jiK) 


n c (cm 3 ) 


N c 


H [31] 


1 


0.0648 [32] 


1.06 x 10" 15 


50 


1.8 x 10 14 


10 9 


Li [33] 


7 


-1.45 [34] 


5.2 x 10~ 13 


0.30 


1.5 x 10 12 


10 3 


Na [35] 


23 


2.8 [36] 


1.90 x 10~ 12 


2.0 


1.5 x 10 14 


10 7 


Rb [37] 


87 


5.4 [38] 


7.3 x 10~ 12 


0.67 


2.2 x 10 14 


2 x 10 6 



Table I. - Comparison of the atomic species that have been Bose condensed. For each cited 
experiment, a is the scattering length, a e \ is the elastic scattering cross section, T c is the tem- 
perature of the onset of BEC, n c is the density of the thermal cloud at the onset, and N c is the 
maximum number of condensate atoms. 



cooling are: 

The elastic collision rate. This determines the rate of thermalization. The low tem- 
perature collision cross section in a Bose gas is determined by the s-wave scattering 
length, a. The elastic collision cross section for identical particles is a e \ — 8na 2 , and 
the elastic collision rate is r c i = ^/2na c \v, where n is the density and v = y/%k&T '/inn. 
Table I shows values of a and a e \ for hydrogen and a number of alkali metal atoms. 
Hydrogen's anomalously small scattering length is conspicuous. In consequence, its elas- 
tic scattering cross section is smaller than that of the alkali metal atoms by a factor of a 
thousand or more, and evaporative cooling proceeds at a correspondingly low rate. 

The cooling path. The ratio of trap depth to thermal energy, rj = Vtrap/ k^T is an 
important parameter in determining the most efficient evaporation path. The rate at 
which atoms with enough energy to escape from the trap are generated is approximately 
proportional to ?7exp(— rj). Each atom carries away energy ~ r/ksT. If rj is large, each 
escaping atoms carries away a great deal of energy, enhancing the efficiency of cooling, 
but because the number of these atoms is small, their evaporation proceeds slowly. If 
■q is too small, the atoms can escape too rapidly for the system to thermalize, and the 
efficiency falls. 

Trap geometry. In a square-well potential, the density of a trapped gas is independent 
of its energy, and the elastic collision rate r c i = \/2nva decreases as v ~ \[T: evaporation 
slows as the temperature falls. In a trap, however, as the cloud cools the atoms are 
compressed into the region of lowest potential, and the density of the gas increases as its 
temperature falls. As a result, in a harmonic trap the evaporation rate increases with 
falling temperature as r e i ~ T~ x . In a quadrupole trap, r c i ~ j'-V 2 . Consequently, 
the functional form of the trap shape is a critical design factor. Because the evaporation 
rate in a trap increases as the temperature falls, one can achieve conditions of runaway 
evaporation in which evaporation, once started, can proceed faster and faster. 

Loss processes. If there were no loss mechanisms in the trapped gas, the time for 
evaporation could be made as long as desired and the magnitude of r e i would be unim- 
portant. However, loss processes are inevitable. In experiments with alkali metal atoms 
the principal loss is scattering of atoms out of the trap by collisions with the warm back- 
ground gas in the cell. The loss rate is independent of the density of the trapped gas. 
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(At very high density, however, loss due to 3-body recombination plays the limiting role.) 
In the cryogenic environment of atomic hydrogen experiments, the major loss process is 
dipolar decay at a rate proportional to density. Because both the dipolar decay rate and 
r cl vary linearly with density, the conditions for runaway evaporation are not achieved. 



6. EVAPORATION TECHNIQUES 

In the initial demonstration of evaporative cooling the atoms were permitted to es- 
cape over a saddle point in the magnetic field at one end of the long cylindrical trap. 
Evaporation is forced by lowering this axial confinement field while simultaneously hold- 
ing the radial confinement fields fixed. Energetic atoms are able to escape out the end of 
the trap. With this method, called saddle point evaporation, it was possible to achieve 
conditions close to BEC in hydrogen, but the cooling power was not adequate to cross 
the barrier. 

The inefficiency of the evaporation process has been explained by Surkov et al. [39] . 
In saddle point evaporation atoms escape only along the z-axis. For an atom to escape, 
it must have a sufficient energy in the axial degree of freedom, E z > Vtrap, where Vt rap 
is the trap depth as set by the saddle point potential. Because only the ^-motion is 
involved, the evaporation is inherently one-dimensional. If the axial and radial motion 
mix rapidly, then all atoms with total energy E > Vt rap can promptly escape. At high 
energy this mixing takes place in our trap because of the coupling of radial and transverse 
co-ordinates. 

As the energy decreases and the trap becomes more harmonic, the mixing time length- 
ens. When it becomes comparable to the collision time, the evaporation rate falls. This is 
because the collision of an energetic atom generally transfers it to lower energy: the atom 
is knocked back into a trapped energy regime. Mixing is governed by the adiabaticity 
parameter, <!> = lu p /lo 2 p , which quantifies the fractional change in the radial oscillation 
frequency per oscillation period as the atom moves along the trap axis (sec cq. 2). For 
$ « 1, the probability of transfering the energy from radial to longitudinal during one 
radial oscillation is ~ $. Using Co p = (dui p /dz)(dz/dt) we obtain 



(7) $ 



0z 



For our trap at the threshold for BEC, $ ~ 10~ 3 . However, the probability for scattering 
during a radial period, 2nT e \/ujp, is about ten times higher. Hence, the evaporation is 
essentially one-dimensional. Surkov et al. have shown that in this situation the cooling 
rate decreases by a factor of about Arj compared to the rate at which evaporation would 
proceed in three dimensions [39]. For hydrogen, such a decrease in the already low 
evaporation rate proves fatal. The effects of one-dimensional evaporation have been 
studied by Pinske et al. [40]. 

The bottleneck of saddle point evaporation is avoided by the technique of radia- 
tive evaporation, "rf evaporation," originally proposed by Pritchard [41], which permits 
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evaporation in three dimensions. A radio-frequency oscillating magnetic field drives 
transitions between the trapped state and some other (untrapped) hyperfine sub-level, 
causing the atom to be ejected from the trap. In hydrogen the trapped hyperfine state 
is (F — l,m — 1), and the rf transition is predominantly to the state (1,0). The reso- 
nance occurs for atoms in a region of the trap where the field magnitude Bq satisfies the 
resonance condition fkv r f — ^b-Bo- Here /Ub-B is the Zeeman splitting between the hy- 
perfine sub-levels. The transition matrix element is TiSIr — [1bB±/\/2, where B± is the 
amplitude of the rf field perpendicular to Bq. The probability of an atom experiencing a 
hyperfine sub-level transition as it traverses the resonance region can be estimated using 
the Landau-Zcner theory. The trapping field is assumed to vary linearly with gradient 
B' in the vicinity of the resonance, and the atom traverses the region at speed v. The 
probability of a two-level atom making a transition as it traverses the resonance region 
is p = 1 — exp(— £) where ( — 2nhilj i / ' \i^,B'v [42]. The hydrogen F = 1 state is actually 
a three level system. Vitanov and Suominen have solved the multilevel Landau-Zcner 
problem [43], and we apply their results to our situation. The probability P m that an 
atom originally in the state (1,1) emerges in the state (l,m) is 



For large ( the atom absorbs two rf photons and emerges in the (1,-1) state, which is 
ejected. In our experiment £ is small, and there is only a small probability of leaving 
the trapped (1,1) state. Transitions are primarily to the (1,0) state. Atoms in this state 
simply fall out of the trap. 

For small (, the probability of being ejected during one radial oscillation is ~ 4(. As 
in the previous discussion on radial and longitudinal mixing, scattering sets a lower limit 
of ( > 10~ 3 . This requires a Rabi frequency of f2# > 2ir x 1 kHz or an rf field strength 
B± > 10~ 7 T. (In other atomic systems dipolar relaxation can impose a more stringent 
requirement on the transition rate [30].) 

To implement rf evaporation in a cryogenic environment methods had to be developed 
to eliminate eddy current losses and rf shielding by the cell. This was accomplished using 
a plastic cell design in which heat transport is provided by a superfiuid helium jacket 
around the cell. The rf field amplitude is typically 7 x 10 T, and fields can be applied 
with frequencies up to 46 MHz. RF evaporation is switched on at a trap depth of about 
1.1 mK, where the sample temperature is typically 120 ^K. 

In addition to driving evaporation, the rf field is used to find the temperature of the 
sample. The temperature is measured by sweeping the rf resonance through the trap 
and measuring the atom ejection rate as a function of frequency [41]. 

The field at the bottom of the trap, the axial bias field 9/hb, is a critical parameter 
because it determines the curvature of the potential minimum and provides the zero point 
for determining the trap depth. If the field is too large the trap is harmonic, the effective 
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volume is large, and the density is small. If the bias field is too low, the condensate 
is strongly confined and suffers a high dipolar decay rate. The optimum bias field is 
approximately ksT/^Q. In such a field the thermal cloud is tightly confined, but the 
condensate can spread out due to interactions. In our experiments the bias field energy 
is about 6 ~ (3/5) A^T 1 . The bias field is measured by applying the rf field at a low 
frequency and sweeping the frequency up until atoms start to leave the trap. 

We have neglected the effect of gravity. For H] at a temperature of 50 /iK the 
gavitational scale height k^T jvng is 4 cm, which is comparable to the vertical size of 
the cloud. For higher mass atoms, which condense at lower temperatures, gravity can 
be important. In such a case the surface of constant B is no longer an equipotential. 
Evaporation occurs primarily from the bottom of the cloud, and becomes one-dimensional 
[30]. 

7. COLD-COLLISION FREQUENCY SHIFT 

At low temperature, in the limit a <C At, only s-wave collisions are important. These 
collisions give rise to a mean interaction energy, and they introduce frequency shifts into 
radiative transitions. Such effects can be analyzed by kinetic theory starting with the 
Boltzmann transport equation [44,45], or described by mean field theory based on the 
pseudopotential [46]. Both approaches give the same result for a homogenous system: 
the mean field energy of an atom in a gas with density n is 

4-7r?i 2 a 

9 Un= nxg 2 (0), 

m 

where the density normalized second order correlation function g 2 (x) is [47] 

(io) fl2 ( x) = _L^ < *| < j (ri _ rj _ x) |*). 

Here ^ is the wave function for the system. (For a Bose gas far from degeneracy, g 2 (0) = 
2.) Because the scattering lengths for 1S-1S and 1S-2S collisions arc not identical, 
the energy to excite an atom to the 25 state from a gas of 15 atoms is shifted by an 
amount [48] 

(f 1) hAv col = (ai 5 _2S - ais_ 13)02(0). 

m 

The frequency shift IS.u co \ is known as a cold collision frequency shift. For a non- 
degenerate gas the two-photon sum frequency is shifted by IS.v co \ = nx, where \ = 
4:h(ais-2S — ais-is)/ m - Once x is known, the density can be determined directly by 
measuring the frequency shift. In addition, a measurement of \ can be used to check the 
theoretical calculations of the scattering lengths. The 15-15 scattering length is known 
accurately from theory: 015-15 = 0.0648 nm [32], and ais-2S has also been computed: 
ais-2S = -2.3 nm [49]. 
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The cold-collision shift of the 1S-2S transition plays a key role in our experiments, 
allowing us to measure the density rapidly in situ. Furthermore, because the density in 
a Bose-Einstein condensate in HI is much higher than the density of the normal gas, the 
cold-collision frequency shift provides an unmistakable signature of the condensate. 

To measure the frequency shift parameter x, a series of line scans were taken at 
different densities as shown in fig. 6a. The initial density was established by monitoring 
the two-body dipolar decay rate, as described above. During successive laser scans of a 
single trapped sample the density decayed because of collisions with helium gas generated 
by heating due to the laser. The area under each photoexcitation curve is proportional 
to the total number of atoms, making it possible to infer the density for each scan. The 
line center for each curve was corrected for the effects of density inhomogeneities due to 
the trapping potential. From a plot of frequency vs. density, like the one shown in fig. 6b, 
the value of x can be determined. From a series of such measurements taken at different 
densities and temperatures, we obtained \ = -3-8 ± 0.8 x 10 10 n Hz cm 3 . The theory of 
cold-collision frequency shifts in an inhomogeneous system is not yet fully understood, 
but assuming cq. 11 is still valid, we deduce ais-2S = — 1.4 ± 0.3 nm, in fair agreement 
with the prediction [50]. 

8. OBSERVATION OF THE CONDENSATE 

Bosc-Einstein condensation in hydrogen was achieved at a temperature of about 50 fiK 
with a density of 1.8x 10 14 cm~ 3 [31]. Its onset was revealed by unmistakable new features 
in the spectrum, shown in fig. 7. In the Doppler-sensitive line there is a narrow peak, 
shifted somewhat to the red of the line center, and there is a similar line to the red of 
the Doppler-free line. (Note: the Doppler sensitive peak was not observed until after 
these lectures.) The shift of these lines reveals a density significantly higher than in the 
normal gas, as expected for the condensate. 

Before BEC was observed in the spectrum there were strong indications that con- 
densation was taking place from a study of the evolution of the density of the thermal 
cloud with decreasing temperature. The density of the non-condensed gas fraction was 
determined from the cold-collision frequency shift; the temperature was inferred from 
the trap depth, set by the frequency of the rf signal. A plot is shown in fig. 8. In these 
observations the value of r\ typically varied between 7 and 5 with decreasing temperature. 
The solid line in fig. 8 is the BEC boundary for i] — 6. Because the normal gas cannot 
exist to the left of the boundary, the density simply falls along the boundary. Once the 
system reaches the boundary, as the temperature is further reduced the normal atoms 
are forced into the condensate. Because the observations are at a laser frequency tuned 
so that only atoms from the normal gas are excited, the condensed atoms, which are 
at such a high density that they are frequency-shifted out of range, are not observed. 
Consequently, the density of the thermal cloud merely tracks the BEC line. 

Because of the condensate's high density, its dipolar decay rate is so high that, in 
isolation, it would disappear in about one second. However, because the condensate is 
continuously fed by the normal gas, its lifetime is several seconds. Its time evolution is 
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Fig. 6. - (a) Series of spectra of a single sample used for a measurement of \. The first scan is 
at the maximum density and exhibits the largest red shift. Subsequent scans, at lower densities, 
are smaller and less shifted. Five of the forty spectra are shown, (b) Shift of the spectrum 
as a function of sample density. The density is inferred from the integrated signal. For this 
particular sample \ = -3.30 ± 0.6 x 10~ 10 Hz cm 3 . From [50]. 



shown in the series of spectral scans shown in fig. 9. As the density in the condensate 
decreases, the red shift decreases and the total signal becomes smaller, finally vanishing. 
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Fig. 7. - Composite 1S-2S two-photon spectrum of trapped hydrogen after condensation, o- 
spectrum of sample without a condensate; •-spectrum emphasizing features due to a condensate. 
The high density in the condensate shifts a portion of the Doppler-free line to the red. The 
condensate's narrow momentum distribution gives rise to a similar feature near the center of 
the Doppler-sensitive line. From [31]. 



9. PROPERTIES OF THE CONDENSATE 

The peak condensate density n o c is found from the red cut-off of the spectrum. As 
shown in fig. 10, n o c = 4.8 ± 1.1 x 10 15 cm~ 3 . The peak mean field energy of an 
atom in the condensate, Un 0tC /k^ = 4irfr 2 an 0tC /kBm = 1.9 jjK, is much larger than 
the energy interval between the radial vibrational states of the trap, Tiiop/k^ = 190 nK. 
Consequently, the shape of the condensate is determined by the mean field energy rather 
than the wave function of the trap's ground state. The density profile, in the Thomas- 
Fermi approximation, is n(r) = n 0jC — V(r)/U. 

With the Thomas-Fermi wavefunction, the total number of atoms in the condensate 
is found by integrating the density over the volume of the condensate. The result is 

2 9/2 7T3/2 5/2 

(12) iV c = Z — ^ %f = (1.1 ± 0.6) x 10 9 atoms. 
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Fig. 8. - Density of non-condensed fraction of the gas as the trap depth is reduced along the 
cooling path. The density is measured by the optical resonance shift, and the trap depth is set 
by the rf frequency. The lines (dash, solid, dot-dash) indicate the BEC phase transition line, 
assuming a sample temperature of (l/5th, l/6th, l/7th) the trap depth. The scatter of the data 
reflects the reproducibility of the laser probe technique and is dominated by alignment of the 
laser beam to the sample. From [31]. 



The oscillation frequencies (eqs. 2) are oj p = 2tt x 3.90±0.11 kHz, and u z = 2n x 10.2 Hz. 
The diameter of the condensate is 15 //m and its length is 5 mm. The huge aspect ratio, 
~ 400, gives the condensate a thread-like shape. 

The fraction of atoms in the condensate is small because the condensate is rapidly 
depleted by dipolar relaxation. The condensate size is limited by the low evaporative 
cooling rate in the normal gas, which supplies cold atoms to the condensate [51]. The 
fraction, / = N c /(N n + N c ), can be found from the integrated area of the normal and 
condensate Doppler-free spectra, taking into account that although the entire condensate 
is in the laser beam, only a portion of the normal gas interacts with the laser. The 
condensate fraction can also be found by comparing the number of condensate atoms, 
determined from n 0tC and the trap geometry, with the number of normal atoms found by 
integrating the Bose occupation function weighted by the trap density of states. (The 
temperature is measured from the Dopplcr-sensitive spectrum of the normal gas.) The 
trap shape cancels in the comparison. The methods are in good agreement at / = 6+ 3 %. 
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Fig. 9. - Time evolution of the condensate. Each spectrum is obtained in 0.67 s. As the 
condensate disappears the spectrum narrows, indicating a reduction in density, and gets weaker, 
indicating a reduction in condensate population. The last panel shows the background count 
rate. From [31]. 



Both methods assume thermal equilibrium, which may not be justified. 

In calculating the density and size of the condensate, we have taken (72(0) = 1 in 
computing the mean field energy (eq. 9), but 32(0) = 2 in calculating n C:Q (cq. 11). 
Although this appears contradictory, it yields a condensate fraction that is consistent 
with the observed intensity ratio of the photoexcitation spectra and with the fraction 
predicted in rcf. [51]. If we consistently take 32(0) = 1, then eq. 12 yields N c = 6 x 10 9 
and / = 25%. If we consistently take g 2 (0) = 2, then N c = 3 x 10 9 and / = 14%. In 
a homogeneous condensate in thermal equilibrium one expects 32(0) = 1, however this 
may not apply under our experimental conditions. The problem clearly requires further 
study. 

The Doppler-free lineshape of the normal gas displays puzzling behavior at the con- 
densation transition. As seen in fig. 11, above the the transition temperature the line 
displays a roughly symmetric shape, as expected for the density distribution in the trap. 
When the condensate is present, the line develops a large asymmetry toward the red. 
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laser detuning [MHz at 243 nm] 

Fig. 10. - Doppler-free spectrum of the condensate (broad feature) and normal gas (narrow 
feature). The dashed line is proportional to the number of condensate atoms at a density 
proportional to the detuning, for an equilibrium density distribution with peak density n 0jC = 
4.8 x 10 15 cm~ 3 in a parabolic trap. From [31]. 



The ratio of normal to condensate volumes is about 10 3 , so that this feature cannot be 
explained simply by penetration of the normal gas into the condensate. 

10. PROSPECTS 

As one expects whenever a novel system is created, the observation of BEC in hy- 
drogen has led to some new questions. In particular, the nature of the photoexcitation 
spectrum from the condensate appears to be more subtle than previously appreciated. 
The cold-collision shift has proven to be an invaluable diagnostic tool, but beyond that it 
has led to an experimental value for an excited-state s-wave scattering length — the first 
such determination to our knowledge. One can conceive of methods for extending such 
measurements to other excited states. 

Perhaps the most dramatic aspect of the hydrogen condensate is its size, more than 
thirty times larger than previous condensates, with prospects for large improvements. 
Thus, hydrogen should be a natural candidate for any application that requires an intense 
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Fig. 11. - Doppler-free spectrum of normal fraction above and below the onset of BEC. The 
symmetric spectrum (above T c , open symbols) suddenly becomes asymmetric (filled symbols) 
when the condensate forms. Temperatures for the three spectra are about 120 /iK (open squares), 
53 piK (open circles), 44 fiK (filled circles). From [31]. 



source of coherent atoms. 

The techniques for achieving BEC in hydrogen evolved over a long time, and the 
apparatus reflects a great deal of history. Vast improvements would be possible if one 
were to start from scratch. In particular, the detection efficiency for Lyman-a photons is 
only about 2 x 10~ 5 , the chief loss in signal being due to the low optical collection solid 
angle, 1.6 x 10~ 2 sr. An apparatus designed for optical access would provide a much 
higher signal rate, permitting a more precise study of the condensate's properties and its 
dynamical behavior. 

The size of the condensate is currently limited by the initial density of hydrogen that 
can be loaded into the trap. Improvements should be possible. On a more speculative 
note, if it were possible to increase the thermalization rate of a hydrogen sample by 
introducing an impurity atom into the trapped hydrogen gas, a major increase in the 
condensate size might be achieved. 
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